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A method for the on the whole minimum-fuel, computer-assisted control of an optionally determined number 
thrusters, arranged as desired on a spacecraft, is to be specified, which method solves this linear optimizati* 
problem within the shortest possible time by reducing the computation expenditure. The method depends or 
simplex algorithm and a suitably provided dual permissible start table. The dual permissible start table is sel 
basis of a previously computed optimum table and the transformation of the actual forces-moments vector b 
to the selected start table. Specific processing steps in the dual simplex algorithm are either not needed at £ 
substituted with processing steps that require considerably less computation expenditure. Finally, a reducec 
used in place of the simplex table. This reduced table consists only of an identification number for a thruster 
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selection. The reduction in computation expenditure is achieved in that specific processing steps are perforr 
ahead of time and the corresponding results are then available in the form of data fields and indices or idenl 
number matrixes in the working memory of the computer. In order to compute the thruster actuation, it is the 
necessary to address and request these results of the initialization part of the method 
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Bes hreibung 

[0001] Die Erfindung bezieht sich auf ein Verfahr n zur insgesamt tr ibstoffminimalen, rechn rgestiitzten Ansteu - 
rung einer beliebig festgelegten Anzahl und beliebig festgelegten Anordnung von Dusen an einem Raumfahrzeug. Die 

5 Dusen konnen dabei diskret Oder auch kontinulerlich ansteuerbar sein. 

[0002] Fur die Ausubung von auBeren Kraften und Momenten an einem , Raumfahrzeug ist der Einsatz von Dusen 
in vieien Situationen notwendig. Da der Einsatz der Dusen Treibstoff verbraucht und dleser unmittelbar EinfluB auf das 
Startgewicht und die Lebensdauer von Raumfahrzeugen und insbesondere von Satelliten hat, kommt der Dusenan- 
steuerung mit in der Summe minimal em Treibstoffverbrauch eine wichtige Bedeutung zu. DOsenanordnungen, die 

10 prinzipiell einen niedrigen Treibstoffverbrauch ermdglichen, eine hdhere Oder gleiche Redundanz bei kleinerer Du- 
sen anzahl aufweisen und/oder durch die Nutzlast am Raumfahrzeug beeinfluBt sind, sind oft unregelm&Bige und 
schiefwinklige DOsenanordnungen. Ausdiesen GrOnden ist die Forderung nach der treibstoffminimalen DGsenansteue- 
rung, insbesondere auch fur beiiebige (schiefwinklige) feste DOsenanordnungen an Raumfahrzeugen von besonderer 
Relevanz. 

15 [0003] Es ist bekannt, daB sich das Problem der treibstoffminimalen Dusenansteuerung auf ein lineares Optimfe- 
rungsproblem zuruckfuhren laBt, fur das der Simplexalgorithmus ein geeignetes und effizientes Ldsungsverfahren 
darstellt (siehe zum Beisplel PJ. Wiktor "Minimum Control Authority Plot: A Tool for Designing Thruster Systems" In: 
Journal of Guidance, Control, and Dynamics, Vol. 17, No. 5, Sept. -Oct. 1994, pp. 4-5). Der Einsatz des Simplexalgo- 
rithmus scheitert jedoch haufig an der benottgten Rechenzeit. So konnen oft die Echtzeitanforderungen zum Beisplel 

20 fur den Einsatz des Simplexalgorithmus innerhalb des Regeltaktes des Lage- und Positionsregelkreises eines Raum- 
fahrzeuges nicht erfullt werden, wei! die besonders in hohen Orbits einsetzbaren Rechner den Simplexalgorithmus 
nicht ausreichend.schnell abarbeiten kdnnen. 

[0004] Eine Vielzah! der sich heute im Einsatz befindenden Raumfahrzeuge (Satelliten) besitzen eine Dusenanord- 
nung, bei der jede DQse zum kdrperfesten Koordinatensystem des Raumfahrzeugs zu einer Achse parallel und damit 

25 Z u den zwet anderen Achsen orthogonal angeordnet ist. Diese orthogonale und achsenparallele Anordnung hat zur 
Folge, daB jede DQse nur Krafte in einem Richtungssinn einer Achse nrichtung und Momente nur in einem Drehsinn 
um die beiden anderen Achsen erzeugen kann. Der Vorteil dieser Anordnung 1st, daB die Steuerung eines jeden Frei- 
heitsgrades des Raumfahrzeuges als starrem Kdrper durch den Einsatz von Jeweils nur zwei Dusen mdglich ist, ohne 
andere Freiheitsgrade zu beeinflussen. Die Ansteuerung der Dusen ist fur diese Dusenanordnung sehr einfach. Der 

30 Nachtei! dieser Dusenanordnung besteht darin, daB mtndestens zw6tf Dusen notwendig sind, um alie Freiheitsgrade 
gteichzeitig und beliebig steuem zu kdnnen. MuB der Ausfall von zwei beliebigen Dusen tolerierbarsein, sind 36 Dusen 
erforderiich. 

[0005] Demgegenuber liegt der Vorteil einer geeigneten schiefwinkligen, also nicht orthogonalen und achsenparal- 
lelen Dusenanordnung im wesentlichen in folgenden drei Punkten: 

35 

1 . Bei frei wahlbarer fester Dusenanordnung sind weniger als 12 Dusen erforderiich. 

2. Fur die Tolerierung des Ausfalls von zwei beliebigen Dusen sind emeblich weniger als 36 Dusen erforderiich. 

3. Bei treibstoffmintmaler Ansteuerung der Dusen ist ein geringerer Treibstoffverbrauch mdglich. 

40 [0006] Der Nachtei I dieser Dusenanordnung besteht darin, daB die treibstoffminimale Dusenansteuerung die Ldsung 
eines linearen Optimierungsprobiems erfordert und der Rechenaufwand fflr die Ldsung dieses Problems mit dem Sim- 
plexalgorithmus oft zu hoch ist. Aus diesem Grunde sind andere Verfahren entwickelt worden, die weniger Rechenzeit 
benotigen und durch die eine mehr Oder weniger trelbstoffminimierte, aber keine insgesamt treibstoffminimale Dusen- 
ansteuerung erreicht werden kann. 

45 [0007] Beim sogenannten "Linear Programming 0 wird das Dusenansteuerungsproblem als lineares Optimierungs- 
problem aufgefaBt und mit einschlagigen Verfahren der linearen Optimierung, wie zum Beisplel den Simplexalgorith- 
mus, getost. Der schon erwahnte Nachteil, der den Einsatz dieses Verfahrens In vieien Fallen unmdglich macht, ist 
der zu hohe Rechenaufwand. 

[0008] Beim Table Look-up n -Verfahren werden alle treibstoffminimalen DOsenanordnungen fur Einheitskrafte und 
so -momente bezuglich einer Achse des kdrperfesten Koordtnatensystems des Raumfahrzeugs berechnet (Ldsen des 
linearen Optimierungsprobiems) und in einer Tabelleabgetegt. Sind im Einsatz Krafte und Momente bezuglich mehrerer 
Achsen gleichzeitig auszuuben, werden die entsprechenden Dusenansteuerungen bezuglich Kraft oder Moment und 
fur jede geforderte Achse und dem geforderten Richtungs-oder Drehsinn aus der Tabelle ausgelesen und mit dem 
entsprechenden Betragswert der Kraft- bzw. Momentenkomponente in dieser Richtung multipiiziert. Den Ansteuerwert 
55 fur jede Duse erhalt man durch die Addition der einzelnen Werte fur die jeweilige Duse. Dieses Verfahren ist sehr 
einfach und sehr schnell, fuhrt im allgemeinen aberzu nicht treibstoffminimalen Dusenansteuerungen. 
[0009] Bei verf inert n und erweit rten "Table Look-up"-Verfahren handelt es sich um Verfahren, die ine Look-up- 
Tabelle anlegen, die Eintrage fur zusatzlich eingefuhrte Achsen enthalt n. Durch diese Verfahren kdnnen treibstoff spa- 
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rende Dusenanst uerungen ermitt Itw rd n. Zugleich nimmt ab rdiellb rsichtlichkeit des Verfahrensab und wegen 
der zunehmenden GroB der Tab lien d r Rechenaufwand erheblich zu. 

[0010] Andere Verfahren beruh n auf der Auswertung und Verbess rung des Erg bnisses, das mit dem Tabl L ok* 
up"- Verfahren ermittelt wird. Das Ziel dieser Auswertung besteht darin, besonders treibstoffverschwendende Dusen - 
5 ansteuerungen zu verhindern bzw. geeignet zu korrigieren. Hier soli exemplahsch ein Verfahren, das auf einem Fuzzy- 
Ansatz beruht, genannt werden (vgl. T.Suzuki, K.Yasuda, S.Yoshikawa, K.Yamada, N.Yoshlda "An Application of Fuzzy 
Algorithm to Thruster Control Systems of Spacecraft" in: ISTS 94-C-16, 19th International Symposium on Space Tech- 
nology and Science, Yokohama, Japan, May 15-24, 1994, pp. 1 - 8). 

[001 1] Eine Treibstoffminimierung der Dusenansteuerung insgesamt hat zum Ziel eine solche Ansteuerung der 00- 
10 sen zu ermitteln, die 

1 . die geforderten Kr&fte und Momente am Raumfahrzeug ausQbt, 

2. fur jede Duse einen nichtnegativen Ansteuerwert Hefert, und 

3. die Summe alter Ansteuerwerte der DOsen mlnimlert. 

15 

[0012] Dies stellt - wie erwahnt - ein lineares Optimierungsproblem dar. Da das Losen des Problems in der Regel 
zu zeltaufwendig ist, kann die Ldsung des quadratischen Optimierungsproblems ais Ersatzproblem durchgefuhrt wer- 
den. 

[0013] Die Losung des Ersatzproblems fuhrt zur Ermittlung derjenigen Dusenansteuerung, die 

20 

1 . die geforderten Krafte und Momente am Raumfahrzeug ausubt, 

2. fur jede DQse einen nichtnegativen Ansteuerwert Hefert, und 

3. die Summe aller quadratischen Ansteuerwerte der Dusen minimiert. 

25 [0014] Dieses quadratische Optimierungsproblem kann uber die Efereitstellung der Moore-Penrose-invertierten Sy- 
stemmatrix der Dusenanordnung am Raumfahrzeug und eines Offsetvektors, der mit Hilfe der singularen Wertezerle- 
gung berechnet wlrd (vgl. z.B. B.Noble, I. W. Daniel "Applied linear algebra" Prentice-Hail International, INC. (1988), 
pp. 338 - 350), geldst werden. Die Ldsung des quadratischen Optimierungsproblems kommt dann als Produkt der 
Moor-Penrose-invertierten Systemmatrix mit dem kommandierten Krafte-Momenten-Vektor und der Addition des Off- 

30 setvektors zustande. Die Addition des Offsetvektors bewirkt, daB im Ldsungsvektor nur nicht negative Ansteuerwerte 
auftreten. Die Berechnung des Offsetvektors wird jedoch dann numerisch immer aufwendiger, wenn die Zahl der an- 
zusteuemden Dusen um zwei oder mehr Qber der Dimension des kommandierten Krafte-Momenten-Vektors liegt. Im 
allgemeinen ist die Ldsung nicht treibstoff -minimal, f&llt aber mit der Ldsung des linearen Optimierungs-problems zu- 
sammen, wenn die Zahl der anzusteuemden DOsen nur um eins grdBer ist als die Dimension des kommandierten 

35 Krafte-Momenten-Vektors (verglefche. z.B. DE 195 19 732 A1) . Weiterhln ist nachtellig, daB langere Einschaltzeiten 
uberproportional "bestraft" werden und daher softener in den Ansteuerungen vorkommen. Fur den Betrieb der DOsen 
und in rem Wirkungsgrad ist es jedoch gunstiger, sie mit langeren Einschaltzeiten zu betreiben. 
[0015] Der Nachteil bestehender Verfahren besteht mlthin darin, daB sle entweder nicht das lineare Optimierungs- 
problem losen und damit im allgemeinen nicht die treibstoffminimale Ansteuerung der Dusen bereitstellen oder aber, 

40 wenn bei den Ansteuerverf ahren das lineare Optimierungsproblem geldst wlrd und der Simplexalgorithmus ausgefflhrt 
wird, dieses ublicherweise zu rechenaufwendig ist. 

[0016] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs angegebenen Art bereitzusteflen, mit 
dem das lineare Optimierungsproblem fur die treibstoffminimale Ansteuerung beliebig fester Dusenanordnungen zur 
gleichzeitigen Realisierung von Kr&ften und Momenten an einem Raumfahrzeug in signifikant kurzerer Zeit geldst 

45 werden kann, als es bisher mdglich ist. 

[0017] Diese Aufgabe wird gemaB der Erfindung durch die im Anspruch 1 gekennzeichneten Merkmale geldst. 
[0018] Das Verfahren nach der Erfindung basiert also vortetlhaft auf einem, dem Simplexalgorithmus angelehnten 
Verfahren, das jedoch aufgrund der konsequenten Ausnutzung der Besonderhetten bei der DQsenansteuerung bei 
Raumfahrzeugen und der Aufgabe der Verwendung des ublichen Simplextableaus zu einer Rechenzeitverkurzung von 

so ca. einer GrdBenordnung fuhrt. 

[0019] Vorteilhafte Ausgestaltungen und WefterfQhrungen des Verfahrens nach der Erfindung sind in den dem An- 
spruch 1 nachfolgenden Patentanspruchen gekennzeichnet. 
[0020] Das Verfahren nach der Erfindung basiert auf folgenden Grundlagen: 

55 a) Das bei der DQsenansteuerung zu Idsende lineare Optimierungsproblem kann wie f otgt aufgeschrieben werden: 

Gesucht ist der Ansteuervektor 1^, der die Ansteuerwerte fur alio DQsen enth&tt und folgende Bedingungen er- 
fullen muB: 
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1 . Durch mOssen di gefordert n Krafte und Mom nt am Raumfahrz ug ausgeubt w rdan. Das 
heiBt, di Gleichung . ^ = b^ muB erfQilt werden, wob i A^y, di Systemmatrlx 1st, In der die drei Kraft- 
und die drei Mom nt nkompon nten fur jed inz Ine Dus in in r Spalt dieser Matrix ing trag n sind. 
Die Systemmatrix hat 6 Zeilen wegen der 6 mdglichen Freiheitsgrade in der Bewegung des Raumfahrzeugs, 
5 namlich 3 Krafte in drei zueinander orthogonalen Achsen und entsprechend 3 Momente urn diese Achsen und 

bei D Dusen D Spalten. 

2 * ISopt - [kopt 1 •• koptj koptoF dan * nur nichtnegative Werte enthalten, 

d.h.: 

10 koptj^°.j = 1 (1)D. 



3. die Summe atler Ansteuerwerte muB minimal sein (Zielfunktion), d.h.: 
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b) Ein Simplextableau besteht aus einem Spaltenvektor der Basisvariablenindizes A*, einem Zeilenvektor der 
Nlchtbastsvariablenlndlzes X^, der Koefflzlentenmatrlx A fwt? f dem Zeilenvektor der auf die Nichtbasfsvariablen 
umgerechneten Zielfunktionskoefflzienten und dem Spaltenvektor der Basisvariablenwerte b fah : 
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lm Verfahren nach der Erfindung kommt das Simplextableau nicht in dieser vollen Form zum Einsatz. 

c) Fur die Umsetzung eines betiebigen 6-dimensionalen Kratte-Momenten-Vektors sind maximal 6 Dusen einer 
45 vorhandenen festen Dusenanordnung erforderlich Oder der kommandierte Krafte-Momenten-Vektor 1st mit der 

vorhandenen Dusenanordnung nicht realisierbar. 

d) Fur jede Dusenauswahl I von 6 aus D Dusen existiert ein Basisvariablenindizesvektor Aj, (I). Die Nummem oder 
Indtzes (1,2,...,D) der in der Auswahl I enthattenen Dusen werden im Basisvariablenindizesvektor Aj, (I) zusam- 
mengefaBt und die Ansteuerwerte dieser Dusen sind die Werte der Basisvariabien des Simplextableaus. Die Werte 

so der Basisvariabien sind entsprechend der Reihenfolge der Dusenlndizes in Aj, (I) im Basisvariablenwertevektor 

btnhCX b.bcom) abgelegt. Zu jedem Basisvariablenindizesvektor Xj, (I) gehort eine Basismatrix Bj, (I). Die Basismatrix 
wird gebildet, indem man diejenigen Spalten der Systemmatrix A^, die zu den ausgewahlten Dusen gehdren, zu 
einer neuen Matrix spaltenweise zusammenstellt. Die Reihenfolge der Spalten in der Basismatrix entspricht der 
Reihenfolge der (ndizes der Dusen im dazugehorigen Basisvariablenindizesvektor A^ (I). Die Nichtbasismatrix B^ 

55 (|) wird aus denjenigen Spalten der Systemmatrix zusammengestellt, die zu den nicht in der Auswahl I enthattenen 

Dusen gehdren. Ihre Indizes sind im Nichtbasisvarlablenindizesvektor X nb (l) zusammengefaBt. Die Reihenfolge 
der Spalten in d r Nichtbasismatrix entspricht der R ihenfolge der Indiz s der Dusen im dazugehorigen Nichtba- 
sisvariablenindizesvektor. Die Ansteu rw rt dies r Dusen sind Null. 
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e) lstfGr inen kommandi rt n Kraft -Moment n-Vektorbfo^ das Optimaltableau Aj, ^ ,Xrtb. opt , A^^, 

und bjab .opt gef und n, dann slnd die Ansteuerwerte der Dusen, die durch di Indlzes in X^.opt adressl rt werd n, 
gleich Null und di j nig n Dtis n, di den Basisvariabl n enteprech n, also durch di Indizes in A^^ adressi rt 
werden, enthaiten als Ansteuerwert den Wert aus der entsprechenden Zeile des optimalen Basisvariabl en werte- 
5 vektors b tah npt . Der Basisvariablenindizesvektor Aj, ^ gibt die optimale Dusenauswahl an. 

f) FQr den Fall der gteichzeitigen AusQbung von alien drei Kraftkompo nenten und alien drei M omenten komp onenten 
ist der Einsatz von maximal 6 der vorhandenen Dusen erforderiich. Es gibt genau 



10 6!(D-6)I 

Moglichkeiten (Kombinationen ohne Wiederholung), 6 aus D DQsen auszuwahlen. Werden die Indizes in den bei- 
den Vektoren der Basisvariablen indizes und Nichtbasisvariablen indizes, Aj, und A^,, stets aufsteigend sortiert, so 
gibt es jeweils genau A(C^) dieser Vektoren. FQr eine teste DQsenanordnung gibt es (nur) eine bestimmte Anzahl 
is von Kombinationen von DQsen, die die Dusenauswahl fur eine Losung des linearen Optimierungs-problems sein 

konnen. Diese Anzahi ist im allgemeinen deutlich kletner als die Anzahl A(C*?). 

[0021] Eine wesentiiche Grundlage des Verfahrens besteht darin, daB stets mit sortierten Basis- und Nichtbasisva- 
riablenindizesvektoren gearbeitet wird. Im weiteren sind die Indizes oder Nummem der Dusen in den Basis- und Nicht- 
20 basisvariablenindizesvektoren z.B. stets der GroBe nach sortiert. Auf diese Weise wird es mdglich, daB einer Dusen- 
auswahl genau ein Indizesvektor, genau eine Koeffizientenmatrix A^ = [a^ kJ l, genau ein Koeffizientenvektor der auf 
die Nichtbasisvariablen umgerechneten Zietfunktion und genau eine Baslsmatrix und Nichtbasismatrix B^, 
zugeordnet werden. 

[0022] Eine welters wesentiiche Grundlage des Verfahrens besteht darin, daB mit el nam reduzierten Tableau gear- 
25 beitet wird. Das reduzlerte Tableau besteht lediglich aus der DQsenauswahlidentiflkationsnummer I und dem Basisva- 
riablenwertevektor bfaH. FQr das Verfahren und die Arbeit mit dem reduzierten Tableau ist neben der Bereitstellung des 
Feldes der Koefflzientenmatrizen alter Optlmaltableaus und der Pivotspaltenindizesmatrtx auch die Bereitstellung der 
DOsenauswahlidentifikationsnummemmatrbc, der Zeilenaustauschlndizesmatrix und der Basis indizesmatrix erforder- 
iich. 

30 [0023] Das Verfahren besteht aus zwei Bestandteiten. Der erste, der tnitialislerungsteil, dient zum einen der Bereit- 
stellung des Feldes der Koefflzientenmatrizen aller Optimaltableaus 



£ - */],] - [/*//]• wobei [a^ kj ] f = dW) - B£\l) ' BJLI) tot 



[0024] Die Dusenauswahlidentrfikationsnummer Ifiuft von 1 bis M ( I = 1 (1 )M ), wobei M die Anzahl aller Dusenaus- 
40 wahlmdgiichkeiten ist, die fur einen belleblgen Krafte-M omenten- Vektor und die vortiegende beliebige und un- 
veranderiiche DQsenanordnung optimal sind. Eine Koeffizientenmatrix eines Optimaltableaus iiegt im Sinne dieses 
Verfahrens dann vor, wenn alle Elemente im zur Koeffizientenmatrix gehdrenden Vektor der auf die Nichtbasisvariablen 
umgerechneten Zietfunktionskoeffizienten 

Etab = Rtab (f^tab) nicnt ne 9 ativ s > nd - FQr die Elemente der Koeffizientenmatrix muB daher gelten: 

45 
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[0025] Zum anderen werden im Initialisierungsteil des Verfahrens folgende AdreBmatrizen: die Pivotspaltenindizes- 
matrix P = [p IK ], die Dusenauswahlidentifikationsnummemmatrix U = [u, J, die Zeilenaustauschlndizesmatrix T = [t (K ] 
und die Basis indizesmatrix L = [1 IK ] bereitgestellt, damit auf die im zweiten Teil des Verfahrens benotigten Daten be- 
55 sonders schneil zugegriffen werden kann und bei der Durchfuhrung der Basisaustauschschritte die Arbeit mit dem 
reduzierten Tableau mdglich wird. Durch die Dusenauswahlidentifikationsnummer 1 = 1,2 M wird jeder Dusenaus- 
wahl, die optimal sein kann, und damit auch jeder Koeffizientenmatrix eines Optimaltableaus eine Identifikationsnum- 
mer zugewiesen. 



5 





EP 0 977 687 B1 



10 



15 



20 



25 



30 



35 



40 



45 



50 



[0026] In der Pivotspalt nindizesmatrix P = [p, K ] wird fur jede Zeiie von jeder in F abgelegten Koeffizientenmatrix 
der Index der dazugehorig n Pivotspalte s G {0, 1 , 2, 3 .... (D - 6)} abgelegt. Fur eine Koeffizientenmatrix und ine 
gewahlte Pivotzeile z ist die Pivotspalt s diejenige, fur die gilt: 



[0027] D.h. die zu einer Zeiie einer Koeffizientenmatrix gehdrige Pivotspalte ist diejenige Spalte, in der das Matrix- 
element der Koeffizientenmatrix in dieser Zeiie negativ ist und der Quotient aus dem Element des Vektors der auf die 
Nichtbasisvariablen umgerechneten Zielfunktionskoeffizlenten in der Spalte und dem Betrag des entsprechenden Ma- 
trixelements aus der Koeffizientenmatrix am kleinsten ist. 

[0028] Mit Hilfe der Pivotspaltenlndizesmatrix ist fQr eine DQsenauswahl I, fur die in F eindeutig die dazugeh6rige 
Koeffizientenmatrix A^h CO abgelegt ist, und gewahlter Pivotzeile z die Angabe der entsprechenden Pivotspalte mit s 
(AiabO), z) = p lz aus der Pivotspaltenindizesmatrix moglich, ohne daB dazu die Koeffizientenmatrix A^l) oder der 
Vektor der auf die Nichtbasisvariablen umgerechneten Zielfunktionskoeffizlenten p^fl) bendtigt wird. 
[0029] Existiert zu einer bestimmten Zeiie keine Pivotspalte, so ist das in der Pivotspaltenindizesmatrix durch eine 
Null (p tz = 0) angegeben. 

[0030] In der Dusenauswahlidentifikationsnummemmatrix U = [u (k ] ist fur jede DQsenauswahiidentifikationsnummer 
I = 1 ,2, .... M und jede Zeiie z = 1 ,2,..., 6 jeder, in F abgelegten Koeffizientenmatrix die Dusenauswahlidentifikations- 
nummer derjenigen DQsenauswahl angegeben, die nach dem Basisaustauschschritt vorliegt, wenn vordem Basisaus- 
tauschschritt die DQsenauswahl I vorlag und die Zeiie z als Pivotzeile fur den Basisaustauschschritt gewahlt wurde. 
D.h., wenn zum Zeitpunkt n die DQsenauswahl l[n] vorliegt und ein Basisaustauschschritt mit der Pivotzeile z[nj durch- 
gefOhrt wird, kann die Identifikationsnummer der neu entstehenden DQsenauswahl ohne Verwendung von A^O) mit 
l[n + 1] = u i[n]z(n] direkt aus der Dusenauswahlidentifikationsnummemmatrix ausgetesen werden. Dies wird dadurch 
moglich, daB durch die Kenntnis des aktuellen Tableaus (bzw. seiner Identifikationsnummer I) und der Wahl der Pivot- 
zeile die Pivotspalte implizit festgelegt ist. Zur Ermittlung des Folgetableaus ist also die explizite Kenntnis der Pivots- 
palte s[n] unnotig. 

[0031] In der Zeilenaustauschindizesmatrix T = [t,J ist fQr jede Zeiie z = 1 ,2, 6 von jeder der in F abgelegten 
Koeffizientenmatrizen der Zeilenindex des Pivotelements im (sortierten) reduzierten Tableau angegeben, das nach 
dem Basisaustauschschritt vorliegt, wenn vor dem Basisaustauschschritt die DQsenauswahl I vorlag und die Zeiie z 
als Pivotzeile gewahlt wurde. 

[0032] in der Basisindizesmatrix L = Q,] ist fQr jede Dusenauswahiidentifikationsnummer der Vektor der Indizes der 
in der DQsenauswahl I enthaltenen DQsen, also derzu I gehdrige Bastevariablenindlzesvektor A^l) = 1 angegeben. 
Die DQsen indizes in den Vektoren \j sind der GroBe nach sortiert. 

[0033] Durch den Initialisierungsteil des Verfahrens wird weiterhtn eine Start-DQsenauswahlidentifikationsnummer I 
[0] = Igfauf und die dazugehdrige invertierte Basismatrix 

Sbitoit^ (SbCstait))" 1 bereitgestellt. Aus numerischen Grunden werden alle Werte in F und ^ , die einen sehr kleinen 
Betrag e, (z. B.: e = 10-*) haben, zu Null Null gesetzt 

[0034] Im zweiten Teil des Verfahrens wird zur Realisierung der im. Krafte-Momentenvektor kommandierten Krafte 
und Momente die optimale Dusenansteuerung, d.h. die Ansteuerwerte der einzuschaltenden DQsen ermitteit. Nur die- 
ser Teil des Verfahrens muB den gesteltten Echtzeitanforderungen genugen. Dieser zweite Teil besteht wiederum aus 
zwei nacheinander auszuf uhrenden Teilen. Zuerst wird ein dual zulassiges reduziertes Tableau als Starttableau ermit- 
teit. Im zweiten Teil wird die erforderiiche Anzahl von Baslsaustauschschritten durchgefQhrt. 

[0035] Ein dual zulassiges reduziertes Starttableau wird wie folgt ermitteit: FQr die Bereitstellung eines dual zulas- 
sigen reduzierten Tableaus ist lediglich die Transformation des aktuellen Krafte-Momente- Vektors b kon) auf die fur den 
Start der Basisaustauschschritte festgelegte DQsenauswahl l[0] = 1^ erf orderllch . Dazu dient die schon berettgestellte 
invertierte Basismatrix B^, ^ BJtt und es gilt: 



[0036] Mit der schon bereitgestellten Dusenauswahlidentifikationsnummer fur den Start l[0] = 1^^ und dem neu 
berechneten Vektor b tah [0] = b tah Rtflrt steht ein reduziertes und dual zulassiges Starttableau [l[0] b^ [0]] bereit. 
[0037] Dann erfolgt die Ausfuhrung der Basisaustauschschritte (zyklischer Verfahrensteil): 

Derfolgendb schrieben V rfahr nst ilmuBwi derholt solang ausg fuhrtw rd n, bisdi optimal DQs nauswahl 
mit den optimalen Anst u rw rten g fund nwurd oder festgestellt wird, daB es fur d naktu llvorlieg nden Krafte- 
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Momenten-Vektor und die vorli gend DQs nanordnung k ine Ldsung gibt. 

[0038] Da von einem dual zulassig n Slmplextabl au in Form des r duzierten Starttableaus ausgegang n wird und 
durch di Basisaustauschschritt imm r wied r nur in Dus nauswahl erzeugt w rden kann, d r n Simplextabl au 
dual zulassig 1st, kann bei der Durchfuhrung der Basisaustauschschritte immer wieder auf die im initiaiisierungsteil 
5 des Verfahrens bereitgestellten Daten insbesondere im Feld der Koeff izientenmatrizen aller optimalen und damit gleich- 
zeitig aller dual zulassigen Simplextableaus und in der Pivotspaltenindizesmatrix zuruckgegriffen werden. 

Erster Schritt: 

10 [0039] Finde einen negativen Wert in b^n]. n in eckigen Klammern bezeichnet die Nummer des Iterationsschrittes. 
Enthalt bfejn] keinen negativen Wert, liegt die Ldsung des Hnearen Optimierungsproblems vor, und es kann der Vektor 
der Dusenansteuerwerte Kopt aufgebaut werden. Enthfilt b^ [n] negative Werte, dann wird eine dieser Zeilen ats Pivot- 
zeile z[n] gewahlt. Es ist gunstig, die Zeile des betragsgroBten negativen Elements in b^n] als Pivotzeile zu wahlen. 

is Zweiter Schritt: 

[0040] Mit Kenntnis der Dusenauswahlidentifikationsnummer l[n], und der im ersten Schritt gewahlten Pivotzeile z 
[n] wird aus der Pivotspaltenindizesmatrix P der dazugehdrige Pivotspaltenindex mit s[n] = p , [nlz(n] ausgetesen. Ist in 
P kein Pivotspaltenindex angegeben (p , [nJzfn j = 0), muB eine andere Pivotzeile ausgewahlt werden. 1st zu keiner der 
20 mdglichen Pivotzellen in der Pivotspaltenindizesmatrix eine Pivotspalte indizlert, extstiert keine Ldsung des Hnearen 
Optimierungsproblems. Sind Pivotzeile und Pivotspalte ermitteft, wird zum dritten Schritt Qbergegangen. 
[0041] Im dritten und im vlerten Schritt wind der eigentliche Basisaustauschschritt durch Ermittlung des n&chsten 
reduzierten Simplextableaus ausgefflhrt: 

25 Dritter Schritt: 

[0042] Mit Hiffe der Koeffizienten in der Pivotspalte s[n] von A^tn] wird aus b^n] der neue Basisvariablenwerte- 
vektor b^.neu berechnet. Die bendtigten Etemente des Pivotspaltenvektors a aus A^Jn] konnen direkt aus dem Feld 
der Koeffizientenmatrizen aller Optimaltableaus mit a^ k s[nJ [n] = f ksln]l[n) ausgelesen werden. 
30 Die Elemente des neuen Basisvariablenwertevektors werden wie folgt berechnet: 
Der neue Wert in der Pivotzeile z[n| wird berechnet nach: 



"tab neu zfrn ~ 7 » 
35 'z[n)sin]t[n\ 

wobei f^j , [n] das entsprechend indizierte Pivotelement aus [n] ist. 
[0043] Die neuen Werte in den anderen Zeilen werden berechnet nach: 

b tabneuk = b tabkl n i ' b tabneuz{n\ ' f ks[n\i[n\ 

wobei: {k= 1 (1)6) aber(k* z[n]) 
[0044] Da der neuberechnete Basisvariablenwertevektor b^^ zum fiktiv angenommenen neuen sortierten Basis- 

45 variablenindizesvektor XJn + 1 ] passen muB, ist Im allgemeinen der neuberechnete Basisvariablenwertevektor ent- 
sprechend umzuordnen. Dazu sind im neuberechneten Basisvariablenwertevektor b fah nm ledigltch die Werte zwischen 
der aktuellen Pivotzeile z[n] und der mit Hilfe von l[n] und z[n] aus der Zeilenaustauschlndizesmatrix T = [t, J auszule- 
senden Austauschzeile z teusch [n] = t i[ n ]z(n] zyklisch zu tauschen. Die Austauschzeile Ztauschl n ] ist die Zeile, wohin der 
Wert des neuberechneten Basisvariablenwertevektors plaziert wird, der in b^ neu in der Pivotzeile z[n] steht. Nach 

so diesem Umordnen des neuberechneten Basisvariablenwertevektors liegt der fur den folgenden Basisaustauschschritt 
benotigte Basisvariablenwertevektor b^n + 1 ] vor. 

Vierter Schritt: 

55 [0045] Mit Kenntnis von l[n] und der Pivotzeile z[n] wird mit l[n + 1 ] - U|[„]z[n] die DQsenauswahlidentifikationsnummer 
ermittelt, die zum neuberechneten und nach aufstelgenden Basisvariablenindizes sortierten Basisvariablenwertevektor 
btab [n + 1] gehort. Der Basisaustauschschritt ist damit. vollzogen, das n ue reduzierte Tabl au [l[n + 1] b^ [n + 1 J] 
ii gt vor, und es kann die Losung des linearen Optimierungs-problems mit dem ersten Schritt des zyklischen V rfah- 
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rensteils fortg setztw rden. 

[0046] Bricht das V rfahren mit einem Basisvariablenw it vektor, d r k in negativ n Werte m hr nthalt, ab, so 
liegt der ptimale 6-dimensi nal Basisvariabl nw rtevektorb^ hnr> = b tHh [n^J v r, und es kann mit Hilf v n l[n end | 
aus der Basisindizesmatrix Lderzum aktuellen Basisvariablenwertevektor gehdrende Vektor der Basisvariablen i ndizes 
s mit A^[n end ] = 1 , tnendl ausgelesen werden. Mit Kenntnis der beiden Vektoren b^n^J und A to [n and ] wird der optimale 
D-dlmensionale Ansteuervektor fOr die Dusen k^ auf gebaut. 

[0047] Fur die Losung des linearen Optimierungsproblems fur unterschiedliche Krafte-Momenten- Vektoren sind im 
allgemeinen auch unterschiedlich viete Basisaustauschschritte n end G {0, 1 , .... n,^ erforderlich. Bei Simulations- 
rechnungen mit einer DQsenanordnung von 12 Dusen wurden in sehr seltenen Fallen bis zu n max = 10 Basisaustausch- 

10 schritte durch gefuhrt, obwohl im Durchschnitt nur ca. 3,5 Basisaustauschschritte notig waren. Bei geeigneter Wahl 
des Starttableaus b teb jrtat J sind jedoch nur In ca. 0,5 % aller Fftlle mehr als 5 Austauschschritte notwendig. Im 
allgemeinen bleibt es ohne merkliche Konsequenzen fOr die Anwendung, wenn die Anzahl der maximal auszufQhren- 
den Austauschschritte auf einen Wert begrenzt wird, bei dem nahezu alte Falle vollstandig bearbeitet werden und in 
den Fallen, in denen doch noch negative Werte im Basisvariablenwertevektor stehen, diese zu Null gesetzt werden. 

is Auf diese Weise kann der Worst-Case-Rechenaufwand weiter gesenkt werden . Anwendu ngen dieser Verfah rensmog- 
lichkeit bei der Simulation in geschtossenen Regelkreisen zeigten keine bemerkbare Verschlechterung der Regetgute 
Oder eine Erhdhung des Treibstoffverbrauchs. 

[0048] Enthait der Basisvariablenwertevektor nach der Ausfuhrung der maximal erlaubten Anzahl von Basisau- 
tauschschritten noch negative Werte, ist es prinzipiell auch moglich, daB keine Losung des linearen Optimierungspro- 

20 blems existiert. Dieser Fall kann dann auftreten, wenn eine Dusenanordnung vorfiegt, die nicht beliebige Krafte und 
Momenta aus dem Beretch aller zugelassenenen Krafte-Momenten-Vektoren erzeugen kann. Damit in diesen Fallen 
keine zu groBen Fehlkrfifte und Fehlmomente erzeugt werden, kann in Abhanglgkeit vom Betrag der negtlven El amenta 
im Basisvariablenwertevektor entschleden werden, ob die negativen Elemente zu Null gesetzt werden und der An- 
steuervektor fQr die DQsen aufgebaut wlrd oder generell keine Ansteuerung der DOsen erfolgt. 

25 [0049] Kommen bei der Abarbettung des Simplexalgorithmus entartete Optimaltableaus vor, dann ist es theoretisch 
mdglich, daB der Algorithmus in einen Endloszyklus I&uft und nicht konvergiert. In Simulationsrechnungen mit unter- 
schiedlichen DQsenanordnungen ist dieser Fall allerdings noch nie eingetreten und wlrd deshalb auch nicht expltzlt 
berQcksichtigt. FOr den dualen Simplexalgorithmus kann man durch die gew&hlte Dusenanordnung sicherstellen, daB 
keine dual entarteten Optimaltableaus existteren. FQr den primalen Simplexalgorithmus, also das Auftreten primal 

30 entarteter Optimaltableaus, ist das so nicht mdglich. 

[0050] Soil bei der Ermittlung der treibstoffminimalen Trlebwerksansteuerung der unterschiedliche spezifische Treib- 
stoffverbrauch cj, j=1 (1)D der einzelnen Triebwerke berijcksichtigt werden, liegt die Zlelfunktion 



E c . K *s 
j w1 j opt j 



40 

[0051] Die Berucksichtigung der Gewichtskoeffizienten ist in dem Verfahren fQr den on-line auszufuhrenden Verfah- 
rensteti ohne Mehraufwand moglich. 

[0052] Falft eine Duse aus, Sndert steh die DQsenanordnung im Sinne dieses Verfahrens. Das Verfahren auf die 
neue DQsenanordnung anzuwenden erfordert, daB alle Daten, die im Initialisierungsteil des Verfahrens bereitgestellt 

45 wurden, fQr die neue DQsenanordnung bereitgestellt werden. 

[0053] Da sich die neue DQsenanordnung nur durch den Wegfall einer DQse von der alten DQsenanordnung unter- 
scheidet, gibt es Mogltchkeiten, die neuen Daten durch eine Aktualisierung der alten Daten zu ermitteln. Durch die 
Bereitstellung eines erweiterten Feldes von Koefflzlentenmatrizen, In dem auch dlejenigen Koeffizlentenmatrizen ab- 
gelegt sind, die durch den Ausfall von DQsen Koeffizlentenmatrizen von Optimaltableaus werden, und entsprechend 

so erweiterter Adre3matrizen , kann der Aktualisierungsaufwand stark reduziert werden und fur bestimmte Matrizen ganz 
entf alien. 

[0054] SchlieBitch besteht die Mdglich keit, die Daten fur alle DQsenanordnungen, die durch den Ausfall einer Duse 
zu erwarten sind, vorab bereitzustellen, und be! Ausfall einer bestimmten DQse dann auf die entsprechenden neuen 
Daten umzuschalten. Bis zum Ausfall einer nachsten DQse mussen fur die dann zu erwartenden DQsenanordnungen 
55 die Daten wiederum bereitgestellt werden, urn dann durch Umschaiten auf den nachsten Dusenausfall reagieren zu 
kdnnen. 

[0055] Auch bei einer V rschiebung des Schwerpunkt s des Raumfahrzeuges kommt es zu einer Veranderung der 
Dusenanordnung im Sinne dieses Verfahrens. I matlgem inen fuhrt das zu in rV randerungderSystemmatrix. Dah r 
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rnussen alle von ihr abgeleit t n Daten n ub rettgest lltwerden. 

[0056] Das zuvor beschriebene Verfahren kann zur weiteren Einsparung von Rechenaufwand und Speicherplatz- 
bedarf noch modifiziert werden: 

5 - In voranstehender Beschreibung wurde stets von der Steuerung alier 6 Freiheitsgrade des Raumfahrzeugs aus- 
gegangen. Das Verfahren ist auch auf die Steuerung von 1 bis 5 Freiheits-, graden anwendbar. Es ist ledigtich die 
Stretch ung bzw. die Nichtberucksichtigung der den nicht zu steuemden Freiheitsgraden entsprechenden Zeilen 
in der Systemmatrix und der von ihr abgeleiteten Vektoren und Matrizen erf orderiich . Desgl eich en ist die Steuerung 
einer sich standig andemden Teilmenge alter 6 Freiheitsgrade mdglich. Liegt eine schiefe Dusenanordnung vor, 

10 bteiben jedoch die nicht gesteuerten Freiheitsgrade im aitgemeinen nicht unbeeinfiuBt. 

Mussen am Raumfahrzeug nicht beliebige Krftfte und Moment e ausgeObt werden, ist der Einsatz einer reduzlerten 
Dusenanordnung mdglich. Das vorliegende Verfahren ist auch ftir derartige Anordnungen anwendbar. Es soltte 
nur gewahrleistet sein, daB die kommandierten Krfifte und Mom ante des Kommandos durch die vorliegende DO- 
15 senanordnung realisierbar slnd. Ist das Kommando nicht realisierbar, gibt das Verfahren an, daB keine Losung 

des linearen Optimierungsproblems existfert. 

In Abhangigkeit von der vorliegenden Dusenanordnung kommt es vor, daB von den im Feld der Koeffizientenma- 
trizen aller Optimaltableaus abgelegten Koeffizientenmatrizen bestimmte Spalten nie als Pivotspatte ausgewahlt 

20 werden. Nach Berechnung der Pivotspaitenindizesmatrix P kdnnen die nicht benotigten Spalten aller Koeffizien- 

tenmatrizen einfach ermtttelt werden. Es ist mdglich, ein Feld aller relvanten Pivotspaltenvektoren V zu erstellen, 
in dem nur die Spalten der Koeffizientenmatrizen aller Optimaltableaus abgelegt slnd, die tatsachlich als Pivot- 
spaltenvektoren ben6tigt werden. Dadurch kann eine signifikante Speicherplatzerspamis erreicht werden. Die 
Adressierung der Matrixspalten im Feld alter relevanten Pivotspaltenvektoren erfordert keinen erhdhten Rechen- 

25 aufwand. Jede Spatte hat im Feld alter relevan-ten Pivotspaltenvektoren ihren Spaltenindex. Wind in einer veran- 

derten Pivotspaitenindizesmatrix, der Ptvotspalteniden-tifikationsnummemmatrix 0 = [o )k | nicht der Spaltenindex 
eines Pivotspaltenvektors in einer Koeffizientenmatrix angegeben, sondem gleich die Pivotspaltenidentifikations- 
nummer als Spaltenindex r im Feld aller relevanten Pivotspaltenvektoren V, dann kann im zweiten Schritt des 
zyklischen Verfahrensteils mit r = o llt gtetch der Spaltenindex im Feld aller relevanten Pivot-spattenvektoren V = 

30 t v krl angegeben werden. Die Ermtttlung der Elemente des neuen Basisvaiiabtenwertevektors im dritten Schritt des 

zyklischen Verfahrensteils andert sich dann wie folgt: 

_ b tabnn\ W 
°tab neu z\n\ ~ v * 
35 V ^n\![n] 

wobei v^pj^j das entsprechend indizierte Pivotelement aus A^n] 1st. Die neuen Werte in den anderen Zeilen 
werden berechnet nach: 

b tab neuk = D tab k t n l * b tabneuz[rr] * V kt[n\ 

wobei: {k= 1(1)6) aber\k* z[nj) Die weiteren Verfahrensbestandteile bleiben unverfindert 

45 - Liegt ein Feld alter relvanten Pivotspaltenvektoren V vor, kann insbesondere bei Dusenanordnungen, die bestimm- 
te Symmetrien aufweisen, festgestellt werden, daB reievante Pivotspaltenvektoren unterschiedlicher Optimalta- 
bleaus numerisch identisch sind. Da das mehrmalige Ablegen gleicher Pivotspaltenvektoren nicht notwendig ist, 
kann emeut weiterer Speicherptatz eingespart werden. Beim sukzesstven Aufbau des Feldes alier relevanten 
Pivotspaltenvektoren im Initialisierungsteii des Verfahrens kann uberpruft werden, ob ein numerisch identischer 

so Pivotspaltenvektor zum gerade ermittelten Ph/otspaltenvektor schon in V abgelegt wurde. Ist dies der Fall, kann 

das emeute Ablegen dieses relevanten Pivotspaltenvektors entf alien und in der zu V gehdrigen Pivotspatteniden- 
tifikationsnummemmatrix 0 wind die Pivotspaltenindentifikationsnummer (Spaltenindex) des schon im Feld alter 
relevanten Pivotspaltenvektoren abgelegten Pivotspaltenvektors eingetragen. 

55 - Besteht das Erfordernis, den Speicherbedarf des Verfahrens noch weiter zu reduzleren, kann eine Verfahrensva- 
riante eingesetzt werden, bei der es nicht erforderlich ist, alle relevanten Pivotspaltenvektoren, sondem nur noch 
Reprasen- , tant nv ktoren im Rechnersp ich rabzul gen.Repras ntantenvektoren sind rel vant Piv tspalten- 
v ktoren, ausdenen mit nur mlnimalem Rechenaufwand weit rer levante Pivotspaltenvektoren ber chnetw rden 
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kdnn n. 

[0057] Jeder Pivotspaltenv ktor a = [a^a 2 000 aj ~° aj 7 gehort zu iner bestimmten Basis, z.B. 

MM = [Vl*1.2~\,.~Vnl T 

und zu einer Nichtbasisvariablen, z.B. = X^. Wenn a ( verschieden von Null ist, laBt sich aus 

10 T 

a(l 1 ,X r )s[a 1 a a ~a,«»a n ] 

derjenige Pivotspaftenvektor 

a nau = ?('2. K) = [aneu.1«imu.2 <**> a rau| oco a rau n ] T 

mit 

20 

a nsu.i = jund a^ik * /) = - a k . a^j 

berechnen, der zur Basis 

^('2)=[^.1 *1.2~*rJ~Vn] T 

und zur Nichtbasisvariablen X nb - X* gehort. So kann aus a(l 1( X r ) mit der Angabe des Austauschindexes i der neue 

30 Pivotspaltenvektor a(l 2 , XJ berechnet werden. 

[0058] Auf diese Weise konnen aus einem z.B. sechsdimensionalen Pivotspaltenvektor bis zu sechs weitere Pivot - 
spaltenvektoren berechnet werden. FOr die Anwendung Im Verfahren zurtreibstoffmlnimalen Dusenansteuerung kann 
man sich jetzt auf das Speichem derjenigen relevanten Pivotspaltenvektoren als Reprasentantenvektoren beschran- 
ken, die erforderlich sind, um alle relevanten Pivotspaltenvektoren wie angegeben zu berechnen. 

35 [0059] Fur die Ermittlung einer mdglichst kfelnen Anzahl von Reprasentantenvektoren lassen sich sehr gunstig Ge- 
netische Algorithmen als Optiomierungsverfahren einsetzen. 

[0060] FQr elne geeignete Adressierung der Reprasentantenvektoren in der Matrix aller Reprasentanten sind statt 
der Pivotspaltenindizes - bzw. Pivotspaltenindentifikationsnummemmatrix einmal eine Reprasentantenidentifikations- 
nummemmatrix R undwelterhin eine AustauschlndizesmatrixS erforderlich. tn der Reprasentantenidentifikationsnum- 

40 memmatrix R wird fur jede DQsenauswahlldenttf Ikationsnummer und Jede Pivotzeiie die Identlfikationsnummer des 
Reprasentantenvektors angegeben, der fur die Berechnung des tatsachlich bendtigten Pivotspattenvektors verwendet 
werden kann. In der Austauschindizesmatrix S wird fur jede DOsenauswahlidentifikationsnummer und jede Pivotzeiie 
der Zeilenindex des Austauschelements angegeben, um aus dem Reprasentantenvektor den tatsachlich bendtigten 
Pivotspaltenvektor berechnen zu kdnn en. Steht an einer Stelle in der Austauschlndexmatrix elne Null, bedeutet dies, 

45 daB der adressierte Reprasentantenvektor unmlttelbar der erforderliche Pivotspaltenvektor ist u nd daher die Umrech- 
nung in den aktuell bendtigten Pivotspaltenvektor entfallt. 

[0061] Da stets mit nach Basisvariablenindizes sortierten Pivotspaltenvektoren gearbeitet wird, muB der zunfichst 
aus dem Reprasentantenvektor berechnete Pivotspaltenvektor noch umsortiert werden (z.B. nach aufsteigenden Ba- 
sisvariablenindizes), um benutzt werden zu konnen. Wird die Sortierung der GrdBe der Basisvariablen vorgenommen, 
50 gibt es bei einem sechselementigen Vektor genau 30 Mdgiichkeiten der Umsortierung, um die erfor-derliche Ordnung 
wieder herzustellen. Bei einer erforderlich en Umsortierung wird stets ein Vektorelement an einer Position im Vektor 
entnommen und auf eine andere Position geschrieben. Die Vektorelemente zwischen den beiden Positionen werden 
entspre-chend verschoben. 

[0062] In einer Sortieridentifikationsnummemmatrix G = [g , k ] wird furjeden Reprasentanten und jede Austauschzeile 
55 die Sortieridentifikationsnummer angegeben, nach der der Vektor umgeordnet werden muB, um den bendtigten Pivot- 
spaltenvektor zu erhalten. 

[0063] Aus einem R prfis ntantenvektor kdnn n durch Baslsaustausch und/od rob n b schrieb n Umordnung 
bis zu 7 x 31 = 217 Vektor n geblldet w rden. Werden alle diese Mdglichkeiten berucksichtigt, kann s aufgrund von 
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Symmetrien in d rTriebw rksanordnung vorkomm n.daBwenig rRepras ntant nv kloren rford rlich sind, als wur- 
dennurdi siebenaus inem Repras ntant nvektorgebildeten r levanten Pivotspaltenvektoren in Betracht gez gen. 
[0064] Eine weitere V rf ahrensvariant rgibt sich daraus, daB auf den Basisaustausch am R prasentantenv ktor 
verzlchtet wlrd und nur die 30 Umordnungsmdgllchkeiten des Reprasentantenvektors berucksichtigt werden. Auf diese 
Weise reprasentiert ein Reprasentantenvektor 31 Vektoren. Aufgrund von Symmetrien in der Triebwerksanordnung 
konnen auch auf diese Weise mehrere relevante Pivotspaltenvektoren generiert werden. Im Vergleich zur tetzten Ver- 
f ahrensvariante werden im atlgemeinen jedoch mehr Reprasentantenvektoren benStigt , jedoch entfallt gleichzeitig der 
Rechenaufwand fur den Basisaustausch am Reprasentantenvektor, und es ist die Austauschindizesmatrix S nicht 
erforderiich. 



Patentansprtiche 

1. Verfahren zur insgesamt treibstoffminimalen rechnergestutzten Ansteuerung einer beliebig festgelegten Anzahl 
und beliebig festgelegten Anordnung von Dusen an einem Raumfahrzeug, 

dadurch gekennzeiehnet, 

daB ein treibstoffminimaler Ansteuervektor (k^) mit ausschlieBlich nicht negativen Ansteuerwerten fur alle 
vorgegebenen Dusen mit Hilfe des dualen Simplexalgorithmus gebildet wird, 

daB dazu jeder Dusenauswahl aus der Gesamtzahl von Dusen am Raumfahrzeug, die die optimal e Dusen- 
auswahl fur die treibstoffminimale Realisierung eines beliebigen Kommandovektors (b^) darstelit, eine Du- 
senauswahlidentifikationsnummer (I) zugeordnet wird, 

daB fur die Ermittlung der optimalen Ansteuerung zunachst der Kommandovektor (bj^) der die auszuubende 
Kraft und/oder Momentenkomponenten enthalt, durch Multiplication mft einer Transformationsmatrix ((Bg^)* 1 ) 
in elnen Startvektor GHWt) transformlert wird, 

daB danach der duale Simplexalgorithmus durchgefuhrt wird, bis der Vektor der Ansteuerwerte (b) nur nicht- 
negative Ansteuerwerte enthalt oder aber die Nichtextstenz einer Losung angezeigt wlrd, 
und daB dabei der duale Simplexalgorithmus realisiert wird, Indem vorab alle Daten, die fur die Abarbeitung 
des dualen Simplexalgorithmus notwendig sind und sich off-line (vorab) berechnen lassen, ermittelt werden, 
in Datenfeldern im Rechnerspeicher abgelegt werden und belm Einsatz des dualen Simplexalgorithmus wah- 
rend des on-line-Betriebes durch geeignete Adressierungen aus dem Rechnerspeicher abgerufen und in der 
Rechnung verwendet werden, wobei nur mit einem reduzierten Simplextableau gearbeitetwird, indem lediglich 
die DQsenauswahlidentifikationsnummer und der dazugehdrige Vektor der Ansteuerwerte (b) gefuhrt werden. 

2. Verfahren nach Anspruch 1 , 

dadurch gekennzeiehnet, 

daB die Transformationsmatrix ((Bstart)) -1 ) die invertierte Basismatrix derjenigen Dusen darstelit, die zur Startdu- 
senauswahlgehdren, diese Dusenauswahl f Or elnen beliebigen Kommandovektor optimal ist, diese Dusenauswahl 
eine Start-Dusenauswahlidentifikationsnummer (1^^ zugeordnet ist und der Startvektor (^b^^) die Ansteuer- 
werte dieser Dusen enthalt. 

3. Verfahren nach einem der Anspruche 1 Oder 2, 

dadurch gekennzeiehnet, 

daB jede Diise am Raumfahrzeug eine Identif ikationsnummer erhalt, und die Werte im Vektor der Ansteuer- 
werte fur eine bestimmte Dusenauswahl stets in einer festgelegten Ordnung anhand der Identifikationsnum- 
mem der Dusen sortiert werden, 

daB bei der Abarbeitung des dualen Simplexalgorithmus nach der Ermittlung einer Prvotzeitennummer und 
mit Kenntnis der Dusenauswahlidentifikationsnummer der aktueil voriiegenden Dusenauswahl die entspre- 
chende Pivotspaltennummer aus einer vorab bereitgestellten Pivotspaltenindizesmatrix ausgelesen wird, 
daB das Nichtvorhandensein einer Pivotspaltennummer zu der voriiegenden Dusenauswahl identif ikations- 
nummer und zu der gewahlten Pivotzeile durch einen besonderen Eintrag in der Pivotspaltenindizesmatrix 
eindeutig angezeigt wird, 

und daB die Nichtexistenz einer Losung dadurch angezeigt wird, daB zu keiner moglichen Pivotzeile in der 
Pivotzeilenindizesmatrix eine Pivotspaltennummer angegeben ist. 

4. Verfahren nach einem der Anspruche 1 bis 3, 

dadurch gekennzeiehnet, 
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daG der Pivotspalt nvektor der Koeffizientenmatrix des aktuellen Simplextableaus, dessen Werte fur di Be- 
rechnung neuer Ansteuerwerte fur die Dusen iner neuen Dusenauswahl bendtigt werden, mit Kenntnis d r 
Dusenauswahl vor dem Basisaustauschschritt und der gewahlten Pivotspaltennummer aus einem v rab b - 
rertgestellten Feld aller moglichen optimalen Koeffizientenmatrizen ausgeiesen wind, 
5 - daG der Vektor der Ansteuerwerte fur die Dusen der Dusenauswahl, die nach dem Basisaustauschschritt 

vorliegt, berechnet wird und die Werte in diesem Vektor entsprechend der Zuordnung zu den Dusen und 
entsprechend der festgelegten Ordnung der Identifikationsnummern fur die Dusen sortiert werden. 

5. Verfahren nach einem der Anspruche 1 bis 4, 
io dadurch gekennzelchnet, 

daG mit Kenntnis der DOsenauswahlidentif ikationsnummer vor dem Basisaustauschschritt und der gewahlten 
Pivotzeiiennummer die DOsenauswahlidentif ikationsnummer derjenigen DQsenauswahl, die nach dem Ba- 
sisaustauschschritt vorliegt, aus einer DCksenauswahlidentifikationsnummemmatrix ausgeiesen wird 
15 und daG nach Voiiiegen eines Vektors der Ansteuerwerte, der nur noch nicht negative Werte enthalt, mit 

Kenntnis der dazugehorigen Dusenauswahl identifikationsnurnmer der Vektor der Dusenindizes aus einer Ba- 
sisvariablenindizesmatrix ausgeiesen wird und die mit ihrer Diisen identifikationsnurnmer indizierten Dusen 
den entsprechenden Wert aus dem Vektor der Ansteuerwerte (b) erhalten sowie die ubrigen Dusen den An- 
steuerwert Null zugewiesen bekommen und dam it der Ansteuervektor (j^pt) bestimmt ist. 

20 

6. Verfahren nach einem der Anspruche 1 bis 5, 

dadurch gekennzelchnet, 

daG nur diejenigen Spalten von Koeffizientenmatrizen dual zulassiger Simplextableaus in einer Matrix aller rele- 
vanten Pivotspaltenvektoren abgelegt werden, die bei der Abarbeitung des Verfahrens tatsachlich als Pivotspal- 
25 tenvektoren ausgewahlt werden kdnnen, 

und daG identische relevante Pivotspaltenvektoren nur einmal in dieser Matrix abgelegt werden. 

7. Verfahren nach Anspruch 6, 

dadurch gekennzelchnet, 

30 

daG fur die Adressierung der relevanten Pivotspaltenvektoren jedem Pivotspaltenvektor in der Matrix aller 
relevanten Pivotspaltenvektoren eine Pi votspalten Identifikationsnurnmer zugewiesen wird und diese in einer 
Pivotspaltenidentifikationsnummemmatrix derart eingetragen werden, daG be! Kenntnis der aktuellen DOsen- 
auswahlidentif ikationsnummer und der aktuell gewahlten Pivotzeiiennummer die Prvotspaltenidentifikations- 
35 nummer des aktuell benStigten Pivotspalten vektors auslesbar 1st, 

daG das Nichtvorhandensein einer Pivotspalte zu einer gewahften Pivotzeile durch einen besonderen Eintrag 
in der Pivotspaltenidentifikationsnummem matrix eindeutig angezeigt wird, 

und daG die Nichtexistenz einer Ldsung dadurch bestimmt wird, daG zu kelner mdglichen Pivotzeile in der 
Pivotspalten identifikationsnummernmatrix eine Pivotspaltenidentifikationsnummer angegeben ist. 

40 

8. Verfahren nach einem der Anspruche 3 bis 7, 

dadurch gekennzelchnet, 

daG statt der relevanten Pivotspaltenvektoren nur bestimmte Pivotspaltenvektoren als Reprasentantenvektoren 
im Rechnerspeicher abgelegt werden, wobei die Reprasentantenvektoren eine moglichst kleine Anzahl von Pivot- 

45 spaltenvektoren darstellen, aus denen die tatsachlich bendtigten relevanten Pivotspaltenvektoren ermittelt werden, 

und daG jeder Reprasentantenvektor eine Reprasentantenidentifikationsnummer zugewiesen bekommt und zur 
Adressierung der Reprasentantenvektoren eine Rep rasentantenidentif ikationsnummern matrix verwendet wird, in 
der fur jede Dusenauswahlidentifikationsnummer und jede Pivotzeiiennummer die Reprasentantenidentifikations- 
nummer desjenigen Reprasentantenvektors eingetragen ist, der fur die Berechnung des erforderlichen Pivotspal- 

50 tenvektors verwendet wird. 

9. Verfahren nach Anspruch 8, 

dadurch gekennzelchnet, 

daG in einer Austauschmatrix fur jede DOsenauswahlidentif ikationsnummer und jede Pivotzeiiennummer der Zei- 
ss lenindex der auszutauschenden Basisvariablen als Austauschelement angegeben ist, urn aus dem Reprasentan- 
tenvektor den bendtigten Pivotspaltenvektor zu berechnen. 

10. Verfahren nach ein m der Anspruche 1 bis 9, 
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dadurch gekennzelchnet, 

daB fur di Mdgtichk it d r V rand rung der Dus nanordnung durch zu rwartend Dusenausffille die ntspre- 
ch nden Dat n vorab berechnet und g sp ichert warden und nach F stst lien ines Dusenausf alls automatisch 
durch den Rechner auf Verwendung der entsprechenden neuen Daten umgeschaltet wird. 

5 

1 1 . Verf ahren nach Anspruch 1 0, 

dadurch gekennzelchnet, 

daB die nach einem Dusenausf all erforderlichen neuen Daten durch eine entsprechende Korrektur und Aktuali- 
sierung der bis zum Dusenausfall verwendeten Daten bereitgestellt werden. 

10 

12. Verfahren nach einem der AnsprQche 1 bis 11, 

dadurch gekennzeichnet, 

daB das Verfahren unabhangig von der Dimension des Kommandovektors (b^™) angewendet wird und dabei fur 
die Berechnung der vorab bereitgesteliten Daten die nichtauszuQbenden Komponenten aus der vollstandigen Sy- 
is stemmatrix der Dusenanordnung eliminiert werden. 



Claims 



20 1. Method of overall fuel-minimal computer-assisted drive control of a desired determined number and desired de- 
termined arrangement of nozzles of a space vehicle, characterised in that 

a fuel-minimal drive control vector (k^ is formed with exclusively non -negative drive control values for all 
predetermined nozzles with the assistance of the dual simplex algorithm, 
25 . for that purpose a nozzle selection identification number (I) is associated with each nozzle selection from the 

total number of nozzles of the space vehicle, which selection represents the optimum nozzle selection for the 
fuel-minimal realisation of a desired command vector (b^™), 

for ascertaining the optimum drive control, initially the command vector (b^n,) containing the force and/or 
moment components to be exerted is transformed into a start vector (b =b start ) by multiplication by a transfor- 

30 mation matrix ((Bs^)" 1 ), 

thereafter the dual simplex algorithm is performed until the vector (b) of the drive control values contains only 
non-negative drive control values or, however, the non-existence of a solution is indicated, 
and in that case the dual simplex algorithm is realised in that firstly all data necessary for the processing of 
the dual simplex algorithm and able to be computed off-line beforehand are ascertained, are filed in data fields 

35 in the computer memory and on use of the dual simplex algorithm during on-line operation are called up from 

the computer store by appropriate addressing and employed in the computation, wherein processing is un- 
dertaken only with a reduced simplex tableau in that merely the nozzle selection identification number and the 
associated vector (b) of the drive control values are recorded. 

40 2. Method according to claim 1 , characterised In that the transformation matrix ((B^f)' 1 ) represents the inverted 
basic matrix of those nozzles which belong to the starting nozzle selection, this nozzle selection is optimal for a 
desired command vector, a starting nozzle selection identification number Oo^) is associated with this nozzle 
selection and the start vector (b = b^g^ contains the drive control values of these nozzles. 

45 3. Method according to one of claims 1 and 2, characterised in that 



each nozzle of the space vehicle receives an identification number and the values in the vector of the drive 
control values for a specific nozzle selection are always classified in a fixed order on the basis of the identifi- 
cation numbers of the nozzles, 

so - in the case of processing the dual simplex algorithm after the ascertaining of a pivot line number and with 

knowledge of the nozzle selection identification number of the nozzle selection actually present the corre- 
sponding pivot column number is read out of a pivot column indice matrix prepared beforehand, 
the non-presence of a pivot column number with respect to the nozzle selection identification number present 
and to the selected pivot line is unambiguously indicated by a special entry in the pivot column indice matrix 

55 - and the non-existence of a solution is indicated in the manner that a pivot column number is indicated with 

respect to no possible pivot line in the pivot line indice matrix. 



4. M thod according to n of claims 1 to 3, characteris d in that 
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th pivot column vector of th coeffici nt matrix of the actual simpl x tableau of which the valu sarerequir d 
for th computation of n w drive control values for th nozzles of a new n zzle s lection is read out f a 
previously prepared fi Id of all possibl optimum coeffici nt matrices with knowl dge of the nozzl s lection 
before the basic exchange step and of the selected pivot column number 
5 - and the vector of the drive control values for the nozzles of the nozzle selection present after the basic exchange 

step is computed and the values in this vector are classified in correspondence with the association with the 
nozzles and in correspondence with the fixed association of the identification numbers for the nozzles. 

Method according to one of claims 1 to 4, characterised in that 

the nozzle section identification number of that nozzle section which is present after the basic exchange step 
is read out of a nozzle selection Identification number matrix with knowledge of the nozzle selection identifi- 
cation number before the basic exchange step and of the selected pivot line number 
and after presence of a vector of the drive control values which contain only non -negative values the vector 
of the nozzle indice is read out of a basic variable in dice matrix with knowledge of the associated nozzles 
selection identification number, and the nozzles denoted by their nozzle identification numbers receive the 
corresponding value from the vector (b) of the drive control values and the remaining nozzles are assigned 
the drive control value zero and thus the drive control vector (K^ is determined. 

20 6. Method according to one of claims 1 to 5, characterised in that only those columns of coefficient matrices of a dual 
permissible simplex tableau are filed in a matrix of all relevant column vectors which can be actually selected as 
pivot column vectors in the processing of the method and that identical relevant pivot column vectors are filed only 
once In this matrix. 



25 7. Method according to claim 6, characterised in that 

for the addressing of the relevant pivot column vectors a pivot column identification number is assigned to 
each pivot column vector in the matrix of all relevant pivot column vectors and these are entered in a pivot 
column identification number matrix in such a manner that the pivot column identification number of the actually 
30 required pivot column vector can be read out with knowledge of the actual nozzle section identification number 

and the actually selected pivot line number, 

the non-presence of a pivot column with respect to a selected pivot line is unambiguously indicated by a special 
entry in the pivot column identification number matrix 

and the non-existence of a solution is determined in the manner that a pivot column identification number is 
35 indicated with respect to no possible pivot line in the pivot column identification number matrix. 

8. Method according to one of claims 3 to 7, characterised in that instead of the relevant pivot column vectors only 
specific column vectors are filed in the computer memory as representative vectors, wherein the representative 
vectors represent a smallest possible number of pivot column vectors from which the actually required relevant 

40 pivot column vectors are ascertained and that each representative vector is assigned a representative identification 

number and a representative identification number matrix is used for the addressing of the representative vectors, 
in which matrix for each nozzle selection identification number and each pivot line number the representative 
identification number of that representative vector is entered which is used for the computation of the required 
pivot column vector. 

45 

9. Method according to claim 8, characterised in that the line index of the basic variables to be exchanged is indicated 
as an exchange element in an exchange matrix for each nozzle selection identification number and each pivot line 
number in order to compute the required pivot column vector from the representative vector. 



so 1 o. Method according to one of claims 1 to 9, characterised in that for the possibility of variation in the nozzle arrange- 
ment by anticipated nozzle failures the corresponding data are computed beforehand and stored and after deter- 
mining a nozzle failure there is automatic switching over by the computer to use of the appropriate new data. 

11. Method according to claim 10, characterised in that the new data required after a nozzle failure are prepared by 
55 an appropriate correction and actualisation of the data used up to the nozzle failure. 

12. Method according to one of claims 1 to 11 , characterised that the method is used independently of the dimension 
of the command vector and in that case for th computation of th previously pr pared data the components 
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which are not to be exerted are eliminated from the complete system matrix of th nozzle arrangement. 
Revendlcatlons 

5 

1. Procede pour la commande globale assistee par calculateur, correspondant a une consommation minimale de 
propergol, d'un nombre et d'une disposition fixes a volonte de tuyeres sur un vehicule spatial, 

caracterise en ce 

10 - qu'un vecteur de commande (k^) correspondant a une quantite minimale de propergol et possedant des 

vaieurs de commande exclusivement non negatives pour toutes les tuyeres predetermines est forme a I'alde 
de I'algorithme du simplexe dual, 

qu*a cet effet un numero (I) d'identification de selection de tuydres est associe a chaque selection de tuyeres 
parmi la totalite de tuyeres du vehicule spatial qui represente la selection optimale de tuyeres pour ta realise- 
is tion, correspondant a une consommation minimale de propergol, d'un vecteur quelconque de commande 

(bjcom). 

que pour la determination de la commande optimale, le vecteur de commande (b^,,,,), qui contient la compo- 
sante de force devant dtre appliquee et/ou la composante de moment devant etre appliquee, est transforme 
par multiplication par u ne matrice de transformation ((Ba^a^e)- 1 ) en un vecteur de demarrage (tHiwmarrage)* 
20 . qu'ensuite I'algorithme du simplexe dual est execute jusqu'a ce que le vecteur d'une valeur de commande (b) 

contienne uniquement des vaieurs de commande non negatives ou bien que la non existence d'une solution 
soit indiquee, et 

que Palgorithme du simplexe dual est mis en oeuvre par le fait qu'au prealable toutes les donnees, qui sont 
necessaires pour le traitement de I'algortthme du simplexe dual peuvent dtre calculees off-line (au prealable) ■, 
25 sont determinees, sont memorisees dans des champs de donnees dans la memoire du calculateur et sont 

appelees a partir de la memoire du calculateur lore de ('utilisation de i'algorithme du simplexe dual pendant 
le fonctionnement on-line au moyen d'adressages appropries et sont utilises dans le calcul, le travail etant 
eff ectue uniquement avec un tableau slmplicial redu It, par le fait que seul le numero d'identification de selection 
de tuyeres et le vecteur, qui est associe, des vaieurs de commande (b) sont pilotes. 

30 

2. Procede selon la revendication 1 , caracterise en ce 

que la matrice de transformation ((B^n^^ga)" 1 ) represente la matrfce de base inversee des tuyeres, qui font partie 
de la selection des tuyeres de demarrage, que cette selection de tuyeres est optimale pour un quelconque vecteur 
de commande, qu'un numero ((^manage) d'identification de selection de tuyeres est associe a cette selection de 
35 tuyeres et que le vecteur de demarrage (b=bddmarraga) contient les vaieurs de commande de ces tuyeres. 

3. Procede selon la revendication 1 ou 2, caracterise en ce 

que chaque tuyere du vehicule spatial recoit un numero d'identification, et que les vaieurs contenues dans le 
40 vecteur des vaieurs de commande pour une selection determinee des tuyeres sont en permanence classees 

selon un ordre fixe, en reference aux numeros d'identification des tuy6res, 

que lors de I'execution de I'algorithme du simplexe dual, apres la determination du numero de ligne-pivot avec 
la connalssance du numero d'identification de la selection actuelle des tuyeres, le numero correspondant de 
coionne-pivot est lu a partir d'une matrice d'indices de colonnes pivots foumtes, 
45 - que Tabsence d'un numero de colonne-pivot en rapport avec le numero present d'identification de selection 

de tuyeres et avec ta ligne-pivot selectionnee soit affichee de facon nette au moyen d'une entree particuliere 
dans la matrice d'indices de colonnes-pivots, et 

que ta non existence d'une solution est indiquee par le fait qu'aucun numero de colonne-pivot n'est indique 
pour aucune ligne-pivot possible dans ta matrice d'indices de lignes-pivots. 



50 



4. Procede selon Tune des revendications 1 a 3, caracterise en ce 



que le vecteur de colonne-pivot de la matrice de coefficients du tableau slmplicial actuel, dont les vaieurs sont 
necessaires pour le calcul de nouvelles vaieurs de commande pour les tuyeres d'une nouvelle selection de 
55 tuyeres, est lu, sur la base de la connaissance de la selection des tuyeres avant I'etape d'echange de base 

et avec la connaissance du numero de colonne-pivot selectionne a partir d'un champ, prepare au prealable, 
de toutes les matrices possibles de coeffici nts, 

que le vecteur d s vaieurs d commande pour les tuyeres d la selection d tuyeres, qui est present apres 
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I'etape d'echange d bas , est calcule Gtqu I svaleurscont nues dans ce vecteur sont classees en fonction 
de I'association aux tuyeres et conformement a I'ordre fixe d s numeros ^identification pour I s tuyeres. 

5. Procede seion I'une des revendications 1 a 4, caracterise en ce 

5 

qu'avec ta connaissance du numero d' identification de selection de tuyeres avant I'etape d'echange de base 
et avec la connaissance du numero seiectionne de ligne-pivot, le numero ^identification de la selection de 
tuyeres, qui est present apres I'etape d'echange de base, est lu dans une matrice de numeros d' identification 
de selection de tuyeres, et 

w - qu'apres I'apparition d'un vecteur des valeurs de commande, qui contient encore seulement des valeurs non 

negatives, et avec la connaissance des numeros d'identification de selection de tuyeres, qui y sont associes, 
le vecteur des indices des tuyeres est lu a partir d'une matrice d'indices de variables de base et les tuyeres, 
qui sont indicees avec leurs numeros d' identification de tuyeres, recoivent la valeur correspondante a partir 
du vecteur des valeurs de commande (b), et que la valeur de commande zero est affectee aux autres tuyeres 

15 et que par consequent le vecteur de commande (k^) est determine. 

6. Procede selon Tune des revendications 1 a 5, caracterise en ce 

que seules les colonnes de matrices de coefficients de tableaux simpliciaux admissible* sont memorisees 
20 dans une matrice de tous les vecteurs importants de colonnes-pivots, qui peuvent etre effectivement selec- 

tionnes en tant que vecteurs de colonnes-pivots, lors de la mise en oeuvre du procede, et 
que des vecteurs identiques Importants de colonnes-pivots sont memorisees seulement une fois dans cette 
matrice. 

25 7. Procede selon la revendication 6, caracterise en ce 

que pour I'adressage des vecteurs importants de colonnes-pivots, un numero d' identification de colonne-pivot 
est affecte a chaque vecteur de colonne-pivot de la matrice de tous les vecteurs importants de colonnes-pivots 
et que ce numero est introduit dans une matrice de numeros ^identification de colonnes-pivots de telle sorte 
30 que, moyennant la connaissance du numero actuel ^identification de selection de colonne et du numero 

actuellement seiectionne de ligne-pivot, le numero d' identification de colonne-pivot du vecteur actuellement 
necessaire de colonnes-pivots peut etre lu, 

que I* absence d'une colonne-pivot en rapport avec une ligne-pivot selection n6 est affiche de facon nette au 
moyen d'une entree particuliere dans la matrice de numeros d' identification d'une colonne-pivot, et 
35 . que ("absence d'une solution est determinee par le fait qu'aucun numero ^identification de colonne-pivot n'est 

indique pour aucune ligne possible-pivot dans la matrice de numeros ^identification de colonnes-pivots. 

8. Procede selon I'une des revendications 3 a 7, caracterise en ce 

40 qu'a la lace des vecteurs importants de colonnes-pivots, seuls certains vecteurs de colonnes-pivots sont me- 

morises en tant que vecteurs representees dans la memoire du calculateur, les vecteurs representees re- 
presentant un n ombre aussi faible que possible de vecteurs de colonnes-pivots, a partir d esq u els sont deter- 
mines les vecteurs importants de colonnes-pivots, qui sont effectivement necessaire, et 
qu'un numero d'identification de representant est affecte a chaque vecteur representatif et que pour t'adres- 

45 sage des vecteurs representees on utilise une matrice de numeros d'identification de representees, dans 

laquelte est produit, pour chaque numero d'identification de selection de tuyere et pour chaque num6ro de 
ligne-pivot, le numero d'identification de representant du vecteur representatif, qui est utilise pour le calcu! du 
vecteur necessaire de colonne-pivot. 

so 9. Procede selon la revendication 8, caracterise en ce que 

dans la matrice d'echange, pour chaque numero d'identification de selection de tuyeres et pour chaque numero 
de ligne-pivot, I'indice de ligne de base devant etre echange est indique en tant qu'element d'echange, pour le 
catcul du vecteur necessaire de colonnes-pivots a partir du vecteur representant. 

55 10. Procede selon I'une des revendications 1 a 9, caracterise en ce 

que pour la possibilite de modifier la disposition des tuyeres en raison de pannes de tuyeres auxqueiles il faut 
s'attendr , les donne s correspondent s sont calculees au prealable t mem risees t, apres la determination 
d'un pann d tuyere, une commutation automatique est effectuee par I calcu lat ur, sur {'utilisation des nouvelles 
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donnees correspondant s. 

11. Procede selon la rev ndication 10, caract6rise n ce 

que les nouvelles don noes necessaires apres une panne de tuyere sont preparees au moyen d'une correction et 
5 d'une actualisation correspondantes des donnees utillsables jusqu'a ta panne d'une tuyere. 

12. Procede selon rune des revendications 1 a 11 , caracterise en ce 

que le procede est utilise independamment de la dimension du vecteur de commande (b^™) et que pour le calcul 
des donnees preparees au prealable, les composantes, qui ne doivent pas etre activees, sont eliminees de la 
10 matrice compldte du systeme de I'ensemble des tuyeres. 
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